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elsewhere In order to represent the position of the crankshaft [1] and [11]. Such 
function has a high level of conrelation with the measured in-cylinder pressure or with 
the measured indicated torque since It effectively represents the motored pressure or 
motored torque information. For the crankshaft speed signal, the relevant signal is the 
crankshaft velocity signal fluctuating around its mean value. Therefore, the general 
correlation function for estimating the tn-cylinder pressure or indicated torque can be 
written as a function of the position function angular speed fluctuation 0, and angular 
acceleration 9, as shown below. 

Estimated Value = ^(/i J,6*) (5) 
ESTIMATION OF IN-CYLINDER PRESSURE 
After the in-cyllnder combustion pressure is estimated based on the crankshaft 
speed measurement, the indicated torque is then calculated accordingly based on the 
estimated in-cylinder pressure and the given engine geometry. The estimation model 
function (basis function) may be set to be the following first-order non-linear model as 
shown In Eq. (6) in order to first estimate the in-cylinder pressure. 

The stochastic estimation approach requires building the cross-correlation 
functions between the estimation quantity (in-cylinder pressure) and the measured 
quantities (three basic variables as well as their cross-terms as shown in Eq. (6)). The 
coefficients, ao through a^, can be obtained by minimizing the mean square difference 
between the measured pressure and the estimated pressure as shown in Eq. (7). 
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As described earlier in Eqs, (2) and (3). taking the partial derivatives with respect 
to each of the coefficients and setting the result equal to zero gives the following cross- 
conpelation matrix system to solve. 



> 



</> </^> </,d> <fod> <6d 
<h> <fi> <//f > <fll> <f^&> 







" </>> " 
























<pU>, 



(8) 



In Eq. (8), the various temis in the matrix represent the cross-conrelations among 
the measured basis variables while the right side of the equation represents the cross- 
correlations between the measured in-cylinder pressure and the measured basis 
variables. These non-linear cross-correlations are pre-computed based on all available 
data at a certain engine operating condition, then the five coefficients are computed 
once for all (cycles and cylinders) at that operating point. Once the coefficients as well 
as these conrelation functions are determined and proper processing has been carried 
out, the estimation procedure reduces down to the simple evaluation of a multivariate 
polynomial form based on the measurements. Therefore, during the estimation phase 
the instantaneous value of the five measured basis variables are used to evaluate the 
simple polynomial equation as shown in Eq. (6) for the desired estimation. Therefore, 
the computational requirements can become very minimal in this approach, and the 
estimation can be achieved In real time with a few computational operations. 

Refemng to Fig. (2), Fig. (2) represents each of the prescribed basis variables 
Including the in-cylinder combustion pressure position function U Based on these 
variables, the in-cylinder pressure was estimated using the basis function described in 
Eq. (6) and the cnoss-conrelation described In Eq. (8). Refening to Fig. (3), Fig. (3) 
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Normaiiz&d RM.S, Error = AlUJll L (9) 

Table (4) illustrates this estimation error for each of the estimations and number 
of resolutions accounted in the computation. Note that the values are averages over all 
engine operating conditions. 

Table 4. Nomialized R,M*S. En"ors for Various Cases 



Estimation 
Type 


Number of Resolutions 


360 


60 


36 


Indicated 


Pressure 


2.694 % 


5.063 % 


3.494 % 


Indicated 
Torque 


Individual 
Cylinder 


3.394 % 


5.810 % 


4.313% 


All 
Cylinder 


6.159 % 


7.603 % 


6.814 % 



ESTIMATION OF INDICATED TORQUE 
The Indicated torque Is estimated directly from the crankshaft speed 
measurements, replacing the two steps procedure of first estimating the in-cylinder 
pressure and secondly calculating the indicated torr|ue accordingly. There are two 
different parts of achieving the indicated torque estimation In this approach. The first 
part is to estimate the individual cylinder torque for each cylinder then calculate their 
summations whereas the other part is to directly estimate the summation of individual 
cylinder torque. 

Basis Function Selection - Various basis functions are investigated in order to 
detemiine the best form of the estimation model for the indicated torque estimation in 
real-time. 

Table 5. Various Basis Functions 



Function 
Number 


Basis Function 


1 
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2 




3 




4 




5 




6 




7 





Considering the estimation accuracy, number of terms, equation order, variable 
selection, etc., several different forms of basis functions were investigated using the 
different resolutions (36, 60, and 360) and all engine operating conditions. Table (5) 
describes each of the basis functions selected from many basis functions that were 
examined. 

Note here that the position function U for estimating the indicated torque is 
different from the previous one used for the in-cylinder pressure estimation. It is 
effectively a normalized motored torque, which can be calculated from the given engine 
geometry, for each individual cylinder as well as summation of all cylinders. 

Coefficient Training - After selecting one of the prescribed basis functions in 
Table (5), the polynomial coefficients were obtained by taldng the same procedures, as 
described in Eqs. (7) and (8). Then, the instantaneous value of the measured basis 
variables or their combinations were used to evaluate each of the polynomial equations 
shown In Table (5) to estimated the desired indicated torque. For instance, choosing 
the basis function 3 would result in the following cross-congelation matrix system. 



</> <jj> <ffi> <f^> <^> 
<fe> <ro> <ff^> <r:^> <fM> 



<}^> <f^> <fi^h> <f^> <^^> 

























<T0> 


.Or. 







(10) 
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